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INTRODUCTION 

Computational fluid dynamics has an increasingly important role in the design 
and analysis of aircraft as computer hardware becomes faster and algorithms become 
more efficient. Progress is being made in two directions: more complex and realis- 

tic configurations are being treated and algorithms based on higher approximations 
to the complete Navier-Stokes equations are being developed. The obvious goal is 
the solution of complete aircraft flowfields with the full Navier-Stokes equations. 

The literature indicates that linear panel methods can model detailed, realistic 
aircraft geometries in flow regimes where this approximation is valid. Solutions 
for nonlinear potential equations for flowfields about nearly complete aircraft con- 
figurations have also been available for some time. Recently, Euler methods have 
progressed to the point of computing flowfield solutions for complete aircraft. 

As algorithms incorporating higher approximations to the Navier-Stokes equations 
are developed, computer resource requirements increase rapidly. Generation of suit- 
able grids becomes more difficult and the number of grid points required to resolve 
flow features of interest increases. Coupling greater grid density with the greater 
execution time per point needed to solve more complex equations creates requirements 
for large memory and long run-time. As a result, early work with the Euler and 
Navier-Stokes equations dealt with one- and two-dimensional flows or three- 
dimensional flow about simple geometries. Recently, the development of large vector 
computers has enabled researchers to attempt more complex geometries with Euler and 
Navier-Stokes algorithms . 

This paper describes the results of calculations for transonic flow about a 
typical transport and fighter wing-body configuration using thin- layer Navier-Stokes 
equations, and about helicopter rotor blades using both Euler /Navier-Stokes equa- 
tions. The two codes employed in the calculations have been developed at McDonnell 
Douglas Research Laboratories and are designated as MDTSL30 (thin-layer Navier- 
Stokes fixed wing-body code) and MDROTH (Euler/Navier-Stokes rotary-wing code). 

Both codes have been fully vectorized for optimum performance on a Cray supercom- 
puter. The codes have also been microtasked on a four-processor Cray X-MP/48 for 
reducing the wall-clock time by judicious use of various Cray software techniques. 

‘[’This research was conducted under the McDonnell Douglas Independent Research and 
Development program. 
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The calculations presented in this paper have been performed on a Cray X-MP/48 
Cray X-MP/14, Cray 2, or Convex C-2. Some of the results are presented in color- 
graphics form produced at a Silicon Graphics IRIS 2400T work station. 

COMPUTATION OF WING/ BODY FLOWFIELDS 


Transonic flowfields about wing/body configurations are calculated by use of 
both the Euler and the Navier-Stokes equations. A code designated MDTSL30 has been 
developed at McDonnell Douglas Research Laboratories to compute flowfields by solv- 
ing the Euler, thin-layer, or slender-layer approximations to the Navier-Stokes 
equations . 

The thin-layer approximation retains viscous terms in one direction. For high- 
Reynolds-number turbulent flow, the dominant viscous effects are the result of dif- 
fusion normal to a body surface. The thin-layer approximation is then suitable for 
geometries where the body can be mapped onto a single plane in the computational 
space, retaining viscous terms normal to this plane. A typical example is an 
isolated wing or fuselage. 

The slender-layer approximation retains viscous terms in two directions, neg- 
lecting only streamwise terms. This approximation is useful for calculating viscous 
flows in regions where two aerodynamic surfaces interact, such as the wing-body 
junction or the wing-tip region. Viscous effects on both the fuselage and wing 
surfaces can be modeled with the slender-layer approximation. 

The MDTSL30 code computes the turbulent flowfield by solving the three- 
dimensional, Reynolds-averaged, thin-layer, or slender-layer approximation to the 
Navier-Stokes equations on body-conforming, curvilinear grids. These equations are 
solved by employing Jameson's finite-volume explicit Runge-Kutta time-stepping 
scheme. Turbulent effects are modeled by employing the Baldwin-Lomax algebraic eddy 
viscosity model. 

The details of the methodology and the test calculations performed to verify the 
code are described by Deese and Agarwal, * Agarwal, Deese, and Underwood, ^ and 
Agarwal and Deese. ^ Microtasking was achieved by employing various Cray software 
tools as described by Booth and Misegades.^ 


Calculations have been performed for three different configurations: (1) the 
0NERA-M6 wing, (2) a typical transport wing-body configuration, and (3) a typical 
fighter wing-body configuration. 

0NERA-M6 Wing : M®=0.8A, 0=3. 0° , Re c =11.7xl0 6 


This standard test case was run on a mesh that had 192 cells wrapped around the 
wing, 32 cells normal to the wing surface and 36 cells in the spanwise direction. 
There were 136 cells on the wing surface in the wrap-around direction and 24 on the 

. y (jj P 

wing surface in the spanwise direction. The value of y = * — at the center of 

the first cell was approximately four with about ten cells in the boundary-layer. 
The grid was constructed by stacking 2-D sectional grids along the wing-span. The 
2-D sectional grids were generated using the MDRL algebraic grid generation code. 
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Figure 1 compares pressure distributions on the wing at five span stations as 
predicted using the Euler and thin-layer models. Agreement with experimental data^ 
is good in the leading edge and trailing edge region. Shock smearing is greater 
than desirable, but can be attributed to the size of the cells in the mid-chord 
region. Cell width is typically 2-3% of the chord in this area with shocks being 
smeared over three or four cells. More grid points in the wraparound direction or a 
better distribution with less clustering in the leading- and trailing-edge regions 
would produce sharper shocks. 


□ Experimental data 

Euler 

Thin-layer 





x/c */ c 

Figure 1. Pressure distributions on an ONERA-M6 wing at various 
spanwise stations; Moo = 0.84, a = 3.0°, Re c = 11.7 x 10 6 , 
192x32x36 mesh. 
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The thin-layer model produces a more accurate prediction of the wing pressure 
distribution than the Euler model. Shock strength is too large and shock position 
is too far downstream when the Euler model is used. The thin-layer model, which 
includes the effects of the boundary layer, captures ' shock location and strength 
more accurately. 

This case was also run on a smaller 140x48x32 grid on a four processor X-MP/48 
to demonstrate the microtasking capability. The grid was limited by the eight 
megaword capacity of the machine. Table I gives a summary of the results. Total 
runtime using the microtasked version of the code was one hour on each of the four 
available processors for a total CPU time of four hours. The residuals dropped 
approximately four orders of magnitude in 2500 iterations at a Courant-Friedrichs- 
Levy (CFL) number of 0.9. The actual speedup in wall-clock time of 3.73 achieved by 
microtasking is very close to the theoretical maximum of 3.77 as determined by 
Amdahl's Law. 


Table I. Performance evaluation of transonic viscous 
wing code MDTSL 30 on Cray X-MP/48. 

• Test Case: ONERA M6 Wing 

M^= 0.84, a= 3.0°, Re c = 11.7 xlO 6 

• Mesh used: 140 (chord wise) x48 (wing-normal)x 32 (span wise) 

• Main Memory Required: 7.5 million words 

• Processing rate on one processor: 2.3 x 10 - * CPU seconds/ 
mesh point iteration 

• 98% of total CPU time spent in four key subroutines - 
EULER, FILTER, TSL1, EMUTURB 

• Theorectical speedup on 4-processors (Amdahl's Law): 

cp — CPU time _ 1 _ ^ ^ 

~ wail clock time” (0.98) f _ o gg) ~ 

4 

• Actual speedup achieved after microtasking = 3.73 

• Processing rate on 4-processors: 6.3 xlO -6 wall clock seconds/mesh 
point/iteration 

• Total CPU time of microtasked code is 5% less than original code 

• Number of iterations N = 2500 

• Root-mean-square residual at N = 1 is 0.106 xlO 4 

• Root-mean-square residual at N = 2500 is 0.4830 

• Reduction/cycle = 0.997 

• Total Cpu time = 14400.00 s 

Transport Wing-Body : ty»=0.76, Of =2° , Re c ~6.4xl0^ 

The grid lines on the surface of a typical transport wing-body configuration are 
shown in Figure 2. A 160(chordwise) x34 (wing-normal) x42(spanwise) mesh for this 
configuration was generated with the 3-D grid-generation code of Chen et al.® There 
are 96 cells on the wing in the chordwise direction and 34 in the spanwise direc- 
tion. The value of y* at the center of the first cell along the wing surface ranges 
from 1.5 to 3 while the value of y + at the center of the first cell along the fuse- 
lage surface varies from 1.0 to 2.5. 
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Figure 2. Typical transport configuration with surface grid lines. 


Figure 3 shows the pressure distribution at five spanwise locations on the wing 
computed with the Euler and thin-layer Navier-Stokes codes and their comparison with 
experimental data. Thin-layer terms over the wing and the fuselage have been in- 
cluded in the calculations. As in the case of the 0NERA-M6, the effect of including 
the thin-layer terms is to weaken shocks on the wing upper surface, move them for- 
ward, and improve the comparison with experimental data. 

A comparison of the slender-layer pressure distributions with thin-layer cal- 
culations and experimental data are shown in Figure 4. As expected there is little 
difference between the pressure distributions over the wing. The major effect of 
including slender-layer terms should be near the wing-fuselage junction. No experi- 
mental data was available for comparison in this region. 

Figure 5 shows pressure contours on the upper surface of the transport wing-body 
as displayed on an IRIS color graphics workstation. The location of the shock on 
the upper surface can be clearly identified. 


Transport Wing-Body : Moo=0.76, a=4.2°, Re c =6.4xl0 6 

Thin-layer calculations were performed for the transport wing-body at an angle- 
of-attack of 4.2°. Navier-Stokes terms were included over the wing but not the 
fuselage to reduce run time, so a 160x34x32 cell grid was used with clustering only 
near the wing surface. The value of y + at the center of the first cell ranges from 
1.5 to 3.0. 

At this angle of attack, a significant region of separated flow appeared on the 
upper surface of the wing. Figure 6 shows the streamwise velocity just above the 
wing surface. The separated flow region runs from about 30% of the local chord to 
the trailing edge starting at about 40% of the span and continuing to near the 80% 
span station. 


Computed pressure distributions are compared with experimental data for five 
spanwise locations in Fig. 7. Near the fuselage and near the wing-tip the predic- 
tions compare well with experimental data. In the separated flow region the code 
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Figure 3. Comparison of Euler and thin-layer Navier-Stokes pressure distributions with experimental data for a typical 
transport wing-body at M oo = 0.76, Re c = 6.4x10*, Of = 2.0°, 160x34x42 mesh. 
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Figure 4. Comparison of slender-layer and thin-layer Navier-Stokes pressure distributions with experimental data for a 
typical transport wing-body at Mq^O. 76, Re c = 6.4x 10 6 , or = 2.0°, 160x34x42 mesh. 
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Figure 5. Pressure distribution, P/P^, on the upper surface of a typical transport wing-body, M oo = 0.76, 

Re c = 6.4 x 10 6 , oc = 2.0°. 

predicts the shock to be too far downstream and stronger than was found experimen- 
tally. This result is consistent with experience using the Baldwin-Lomax turbulence 
models for transonic airfoil calculations under separated flow conditions. 

A color contour map of the wing-body upper surface pressure is shown in Fig. 8. 
Again the shock location is clearly defined. 

Fighter Wing-Body : Moo=0.90, a=4.8°, Re c =5.4xl0^ 

The grid lines on the surface of a typical fighter wing-body configuration are 
shown in Fig. 9. The mesh was generated using the 3-D grid-generation code of Chen 
et al.^ The grid has 144 cells in the chordwise direction, 34 cells normal to the 
wing, and 32 cells in the spanwise direction. Only the thin-layer terms near the 
wing have been included for this configuration so there is no grid clustering near 
the fuselage to accommodate viscous effects. There are 96 cells on the wing surface 
in the chordwise direction and 24 in the spanwise direction. The value of y + at the 
center of first cell along the wing surface ranges from 1.0 to 3.0. 

Figure 10 shows the pressure distribution at five spanwise locations on the wing 
computed with the Euler and thin-layer Navier-Stokes codes and their comparison with 
experimental data. Again the viscous effects tend to weaken the wing upper surface 
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Figure 6. Velocity contours, u/VF^/p^, on the upper surface of a typical transport wing-body at 
Mqo = 0.76, Re c = 6.4 x 10 6 , and a = 4.2°. 

shocks and move them forward, although the effect is less pronounced for this 
configuration than for the transport wing-body discussed earlier. 

Pressure contours on the upper surface of the fighter wing-body are displayed in 
Fig. 11. 


COMPUTATION OF HELICOPTER ROTOR FLOWFIELDS IN HOVER AND FORWARD FLIGHT 

An Euler/Navier-Stokes code, designated MDROTH, has been developed at McDonnell 
Douglas Research Laboratories for calculating the transonic flowfield of a 
multibladed helicopter rotor in hover and forward flight. 

The code solves the three-dimensional Euler or Navier-Stokes equations in a 
rotating coordinate system on body-conforming curvilinear grids around the blades. 
The equations are recast in absolute-f low variables so that the absolute flow in the 
far field is uniform, but the relative flow is non-uniform. Equations are solved 
for the absolute-f low variables by employing Jameson's finite-volume explicit Runge- 
Kutta time-stepping scheme. Rotor wake effects are modeled in the form of a correc- 
tion applied to the geometric angle-of-attack along the blades. This correction is 
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Figure 7. Comparison of thin-layer Navier-Stokes pressure distributions with experimental data for a typical transport 
wing-body at M oo = 0.76, Re c = 6.4xl0 6 , a = 4.2°. 
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Figure 10. Comparison of Euler and thin-layer Navier-Stokes pressure distributions with experimental data for a typical 
fighter aircraft wing-body at Mq^O. 90, Re c = 5.4xl0 6 , a = 4.8°, 144x34x32 mesh. 


ORIGINAL PAGE IS 
OF POOR QUALITY 












Figure 11. Pressure contours, P/P^, on the upper surface of a typical fighter configuration; Mq^O.9, 

Re c = 5.4x 10 6 , a = 4.8°. 

obtained by computing the local induced downwash with a free-wake analysis program. 
The details of the methodology are described by Agarwal and Deese.^“9 

A set of test calculations was performed by Agarwal and Deese^ - ^ to verify the 
code for a model rotor in hover and forward flight at various collective pitch 
angles. The model rotor has two, untwisted, untapered blades of aspect ratio equal 
to six and a NACA 0012 airfoil section. Experiments on this rotor have been per- 
formed by Caradonna and Tung^^ at NASA Ames Research Center for a range of blade tip 
Mach numbers, M t , and collective pitch angles, # c . 

Like MDTSL30, the MDROTH code has been fully vectorized for optimum performance 
on a single-processor Cray X-MP and microtasked on a four-processor Cray X-MP/48 to 
reduce the wall-clock time by judicious use of Cray software techniques (Booth and 
Misegades^). The actual speedup in wall-clock time achieved by microtasking (3.73) 
is very close to the theoretical speedup possible (3.77), as was the case for 
MDTSL30. 

Figure 12 shows the main features of the flowfield of a two-bladed rotor in 
hover, with the imbedded finite-difference grid and the coordinate system. The 
flowfield is characterized by transonic shocks, complex vortical wakes, and blade- 
vortex interactions. 
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Figure 12. Schematic of a helicopter rotor flowfield in hover. 


Hover Flowfield Calculations - Euler Solutions 

Computations have been performed for the flowfield of a model helicopter rotor 
in hover. The rotor has two blades of aspect ratio equal to six which are un- 
twisted, untapered and have the NACA 0012 airfoil section. Experiments on this rotor 
have been performed by Caradonna and Tung^O for a range of blade tip Mach numbers 
and collective pitch angles, 9 C . Euler calculations are presented for a hovering 
rotor at a tip Mach number of M t = 0.44 with a collective pitch angle of 9 C = 8.0° 
and at M t = 0.88 with 9 C - 8.0°. 

Rotor-wake effects are modeled in the form of a correction applied to the 
geometric angle of attack of the blades. The free-wake analysis program HOVER of 
the McDonnell Douglas Helicopter Company was used to calculate the induced downwash 
velocity v^ at each span station. For both the cases presented here, the induced 
angle of attack a ± can be taken as approximately 3.8° over 95% of the span so that 
the effective angle of attack becomes o e = 9 C - = 4.2°. 
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These computations were performed on a 98 (chordwise) x33 (blade-normal) 
x21(spanwise) mesh. A typical computation requires two million words of main 
memory, and 2.1 # 10“^ seconds of CPU time per mesh point for each iteration. A 
converged solution required approximately six minutes of CPU time on a single 
processor of Cray X-MP/48 using the fully vectorized version of MDROTH. 

Comparisons of computed pressure distributions with experimental data (Caradonna 
and TunglO) are shown in Fig. 13 for a hovering rotor at a tip Mach number of 

M t * 0.44 and collective pitch angle 6 C = 8° , and in Fig. 14 for M t = 0.88 and col- 
lective pitch angle 6 C = 8°. Both cases indicate good agreement; further improve- 
ment can be achieved by refining the wake model. 

Hover Flowfield Calculations - Navier-Stokes Solutions 


Navier-Stokes calculations for the model helicopter rotor in hover have been 
performed at M t = 0.44, 0 C » 8.0° and M t = 0.61 and 8 C = 8.0°. A 97x33x33 mesh with 
grid clustering near the blade surface was used to resolve the viscous layer. The 
effective angle of attack is 4.2° for both cases. 

Figure 15 shows pressure distributions for the M t = 0.44 case, while M t = 0.61 
results are presented in Fig. 16. Including the viscous effects had little effect 
on the surface pressure for those two cases. No other data are available for 
comparison. 

Forward Flight Flowfield Calculations - Euler Solutions 

Forward flight calculations have been performed for the model NACA 0012 rotor, 
the 0LS rotor, the MDHC 500E rotor, and the AH-64 Apache rotor. All forward flight 
computations have been performed on a 97x33x33 mesh. A typical calculation for the 
full rotor cycle from f=0° to ^=360° requires 12 to 16 hours of CPU time on a Cray 
X-MP/14; approximately 80 iterations are required for one degree movement in the 
azimuth direction. The code is run in time-accurate mode starting with f-0 and 
freestream conditions. 

Figures 17 and 18 show a comparison of the computed pressure distributions on 
the model helicopter rotor in forward flight with experimental data 1 ^ for a location 
near the blade tip, and 

(a) tip Mach number M t = 0.7, advance ratio /J=0.3, and collective rotor pitch angle 

0 c =o°, 

(d) tip Mach number M t =0.8, advance ratio fi= 0.2, and collective rotor pitch angle 

0 C = 0° , respectively . 

In both cases, the computed results show good agreement with the experimental 
data. 

Computations were performed for the OLS rotor blade at M t = 0.63, [l = 0.30 and a 
collective pitch angle 9 C = 0. Results from the NASA Ames Research Center code 
CAMRAD were used to determine the effective angles-of-attack along the OLS rotor 
blade at various azimuth angles. Pressure distributions computed with the Euler 
code are compared with inflight data at the 95% span station as shown in Fig. 19. 

The calculations show good agreement with the experimental data. 
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Euler calculations (with wake correction) 

° a Experimental data 
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Euler calculations (with wake correction) 

O □ Experimental data 


Figure 14. Surface pressure distributions on a lifting 

rotor in hover; M t = 0.88 , 0 C = 8°, AR = 6.0, 
untwisted, untapered, NACA 0012 blade, 
97x33x21 mesh. 
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Figure 13. Surface pressure distributions on a lifting 
rotor in hover; M t = 0.4, 0 C =8°, AR = 6, 
untwisted, untapered, NACA 0012 blade, 
97x33x21 mesh. 














Figure 15. Surface pressure distribution on a lifting rotor in hover; 
M ( = 0.44, Re c = 2xl0 6 , 0 C = 8°, AR = 6.0, untwisted, 
untapered, NACA-0012 blade, 97x33x33 mesh. 


Figure 20 shows the pressure distribution, using the Euler code, computed on 
the 500-E rotor blade for a location near the blade tip. The calculation was 
performed for M t « 0.58, /l = 0.312, and 9 C = 0° using effective angles-of-attack 
given by CAMRAD. Mach number contours on the blade upper surface at various rotor 
angles are shown in Fig. 21. No experimental data are available for this case. 

The final forward flight results to be presented are for a helicopter rotor. 
Calculations were performed at a tip Mach number, M t = 0.65 and ft - 0.33. Effective 
angles-of-attack were determined using CAMRAD. Pressure distributions at the 95% 
span station are shown in Fig. 22. Mach contours on the blade upper surface are 
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Figure 16. Surface pressure distributions on a lifting rotor in hover; 
M, = 0.612, Re c = 2.67x 10 6 , 0 = 8°, AR = 6, untwisted, 
untapered, NACA-0012 blade, 97x33x33 mesh. 


displayed at various azimuth angles in Fig. 23. An expanded view of the Mach con- 
tours near the tip at f = 90° is shown in Fig. 24. A significant region of super- 
sonic flow is apparent on the upper surface. Experimental data are not available for 
comparison. 

CONCLUDING REMARKS 

An extensive program to develop advanced CFD codes based on the Euler and 
Navier-Stokes equations is being conducted at the McDonnell Douglas Research 
Laboratories for applications to helicopters, transport and fighter aircraft, and 
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Figure 17. Pressure distribution on a model rotor-blade of NACA 0012 airfoil section 
and AR = 6 in forward flight; tip Mach number M t = 0.7, advance ratio 
H = 0.3, collective pitch 0 C = 0, 96x32x32 mesh. 



Figure 18. Pressure distribution on a model rotor-blade of NACA 0012 airfoil section 
and aspect ratio 6 in forward flight; tip Mach number M ( = 0.8, advance 
ratio ft- 0.2, 96x32x32 mesh. 
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Figure 19. Pressure distribution on the upper surface of the OLS rotor-blade in forward flight at 95% span location 
tip Mach number M, = 0.63, advance ratio u = 0.30, collective pitch 0„ = 0, 96 X 32 x 32 mesh. 
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Figure 20. Computed pressure distributions on the upper and lower 

surface of the MDHC 500E rotor-blade in forward flight at 
95% span location; M| = 0.58, /* = 0.312, 96x32x32 cell mesh. 


missiles and hypersonic vehicles. Representative calculations about a transport 
wing-body, a fighter wing-body, and a helicopter rotor clearly demonstrate that the 
state-of-the-art in CFD has progressed to the point that turbulent -flow calculations 
about complete fixed- and rotary-wing aircraft configurations may be achieved in the 
near future. 

Efficient use of large computers, including multiple-processor facilities, is 
receiving special attention. 
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Figure 21. Mach number contours on the upper surface of the 500E rotor-blade; M ( = 0.58, \x = 0.312. 
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Figure 22. Computed pressure distribution on the upper and lower 

surface of a helicopter rotor-blade in forward flight at 95% 
span location; M t = 0.65, fi= 0.33, 96x32x32 cell mesh. 
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Figure 24. Mach number contours on the upper surface of a helicopter rotor-blade; M, = 0.65, /u = 0.33 
and 'I' = 90°. 


ORIGINAL PAGE 

C6i0R photograph 


5 




REFERENCES 


1. Deese J.E. and Agarwal R.K. (1987), Navier-Stokes Calculations of Transonic 

Viscous Flow About Wing-Body Configurations, AIAA paper 87-1200, AIAA 19th 
Fluid Dynamics, Plasma Dynamics, and Laser Conference, Honolulu, Hawaii, June 
1987. 

2. Agarwal R.K. and Deese J.E., and Underwood R.R. (1985), Computation of 

Transonic Viscous Wing-Body Flowfields Using Unsteady Parabolized Navier- 
Stokes Equations, AIAA paper 85-1595, AIAA 17th Fluid Dynamics, Plasma 
Dynamics, and Laser Conference, Cincinnati, Ohio, June 1985. 

3. Agarwal R.K. and Deese J.E. (1984), Computation of Viscous Airfoil, Inlet, and 

Wing Flowfields, AIAA paper 84-1551, AIAA 16th Fluid Dynamics, Plasma 
Dynamics, and Laser Conference, Snowmass, Colorado, June 1984. 

4. Booth M. and Misegades K. (1986), Microtasking: A New Way to Harness 

Multiprocessors, Cray Channels, vol. 8, pp. 27-29. 

5. Schmitt V. and Charpin F. (1973), Pressure Distributions on the 0NERA-M6 Wing 

at Transonic Mach Numbers, AGARD-AR-138, chapter B-l. 

6. Chen L.T., Vassberg J.C., and Peavey C.C. (1984), A Transonic Wing-Body 

Flowfield Calculation with Improved Grid Topology and Shock- Point Operators, 
AIAA paper 84-2157, AIAA 16th Fluid Dynamics, Plasma Dynamics, and Laser 
Conference, Snowmass, Colorado, June 1984. 

7. Agarwal R.K. and Deese J.E. (1987), Euler Calculations for Flowfield of a 

Helicopter Rotor in Hover, J. Aircraft, vol. 24, pp. 231-238. 

8. Agarwal R.K. and Deese J.E. (1987), An Euler Solver for Calculating the 

Flowfield of a Helicopter Rotor in Hover and Forward Flight, AIAA paper 
87-1427, AIAA 19th Fluid Dynamics, Plasma Dynamics, and Laser Conference, 
Honolulu, Hawaii, 8-10 June, 1987. 

9. Agarwal, R. K. and Deese, J. E. (1988), Navier-Stokes Calculations of the 

Flowfield of a Helicopter Rotor in Hover, AIAA paper 88-0106, AIAA 26th 
Aerospace Sciences Meeting, Reno, Nevada, January 1988. 

10. Caradonna F.X. and Tung C. (1981), Experimental and Analytical Studies of a 
Model Helicopter Rotor in Hover, Vertica, vol. 5, pp. 149-161. 


545 



